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1. Introduction 

The nonlinear Schrodinger equation with additional potential f/(x), 

i^t = -A^ + f/(x)^ + a|*P^, (1) 

A = — — — a-±l 
dx"^ dy"^ dz"^ ' 

arises in many physical applications including models of optics |T], plasma physics [2] 
and theory of ultracold gases [3]. In the last context, Eq.(Il]) (called the Gross-Pitaevskii 
equation, GPE) appears as one of the basic equations to describe the phenomenon of 
Bose-Einstein condensation (BEG) in so-called mean-field approximation. In this case 
\E'(t,x) means the macroscopic wave function of the condensate, cr = 1 corresponds to 
the case of to repulsive interparticle interactions and a = —1 - to the case of attractive 
interactions. The function U{x.) describes the potential of the trap to confine the 
condensate. In particular, magnetic trap has been modelled by the parabolic potential 
f/(x) = |xp and optical trap has been described by the potential which is periodic with 
respect to one or several variables [H El [6]. 

An important class of solutions for Eq.([T]) are stationary nonlinear states defined 
by the anzats 

^(t,x) = e-^'^V(x). (2) 

The parameter u in terms of BEG corresponds to the chemical potential. The function 
ip{x.) solves the equation 

Ailj + {io- U{^))iIj - (t|V^| V = 0. (3) 

It is known that Eq.(l3]) describes a great variety of nonlinear objects. In particular, it 
has been found that real ID- version of Eq.Q 

V'xx + itu- U{x))i, - = (4) 

with the model cosine potential 

U{x) = Acos2x (5) 

describes bright and dark gap solitons [3 [3 [9l [10], nonlinear periodic structures 
(nonlinear Bloch waves) [3, E], domain walls [I2], gap waves [13] and so on. Some 
interesting relations between various nonlinear objects described by Eq.(jl]) have been 
observed. In particular, in papers [TH [T5] the composition relation between gap solitons 
and nonlinear Bloch waves was established: it has been observed that a nonlinear Bloch 
wave can be approximated by an infinite chain of narrow gap solitons (called fundamental 
gap solitons, FGS), each localized in one well of the periodic potential. In [16] this 
principle has been applied to the case of more general nonlinearity. It is worth noting 
that the results of [13] also can be interpreted in a similar sense, since the gap waves 
discovered in [13] can be regarded as compositions of finite number of FGS. 

In the present paper, we address the problem of description of nonlinear states 
covered by Eq.(jl]) in the case of repulsive interactions, cr = 1, i.e for the equation 

+ U{x))iP -^3 = 0. (6) 
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We argue that if the periodic potential U{x) satisfies some conditions, then all the 
solutions of Eq.(|6]) defined at the whole M can be put in one-to-one correspondence 
with bi-infinite sequences of integers n = 1, . . . ,N (called codes). The correspondence 
is a homeomorphism for properly introduced topological spaces. Each of the integers 
n is "responsible" for the behavior of the solution ip{x) on one period of the potential 
U{x). From this viewpoint, the solutions ip{x) may be regarded as compositions of FGS 
localized in the wells of the periodic potential and taken with a proper sign. So, the 
coding technique gives a unified approach to describe both gap solitons and nonlinear 
Bloch waves and generalizes (and justifies) the composition relation of [TUITS]. In order 
to conclude that for a given U (x) the coding is possible, one has to verify numerically 
three Hypotheses formulated in Section HI As an example, we applied this method to 
the case of model periodic potential ([5]) and present the regions in the parameter plane 
{u, A) where all nonlinear states can be encoded with bi-infinite sequences of integers. 

Our approach is based on the following observation: the "most part" of the solutions 
for Eq.(l6]) are singular, i.e. they collapse (tend to infinity) at some finite point of 
real axis. The set of initial data at x = for non- collapsing solutions can be found 
numerically by properly organized scanning procedure. Then we study transformations 
of this set under the action of Poincare map using methods of symbolic dynamics. A 
similar idea was used to justify a strategy of "demonstrative computations" of nonhnear 
modes for ID GPE with repulsive interactions and multi-well potential pTl ITS]. This 
allowed to find numerically all the localized modes for Eq.Q with single- well and double- 
well potentials and to guarantee that no other localized modes exist. 

The paper is organized as follows. In Section [2] we introduce some notations and 
definitions which will be used throughout the rest of the text and make some assertions 
about them. In Section [3] we formulate a theorem (Theorem [3ll) which gives a base for 
our method. Section H] contains an application of Theorem [31 1 to the case of Eq.([6]). 
The main outcome of Section H] is formulated in the form of three hypotheses. These 
hypotheses provide sufficient conditions for the homeomorphism mentioned above to 
exist and should be verified numerically. In Section [5] we apply this approach to the 
case of the cosine potential ([5]). Section [6] includes summary and discussion. 

For the sake of clarity all the proofs are removed from the main text to Appendices. 

2. Bounded and singular solutions 

2.1. Some definitions 

In what follows we refer to a solution ip{x) of Eq.([6]) as a singular solution if for some 

X = Xq 

lim iIj{x) = +00 or lim tplx) = — oo. 

x—>xo x—^xo 

In this case we say that the solution ip{x) collapses at xq. Also let us introduce the 
following definitions: 
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Collapsing and non- collapsing points: A point {ipO:i^o) plane = {jp^^j') is 

• L-coUapsing forward point, L > 0, if the solution of Cauchy problem for Eq.(l6]) 
with initial data ipi^) = i^o, i^xi^) = i^'o collapses at value x = Xq and < Xq < L; 

• L-non-coUapsing forward point, L > 0, if the solution of Cauchy problem for Eq.(l6]) 
with initial data ipi^O) = ipQ, V'xlO) = V'o does not collapse at any value x = Xq, 
< xo < L. 

• L-coUapsing backward point if the corresponding solution of Cauchy problem for 
Eq.(l6]) collapses at some value x = —xq and < xq < L; 

• L-non-coUapsing backward point if the corresponding solution of Cauchy problem 
for Eq.(l6]) does not collapse at any value x = —xq, < xq < L; 

• oo-non-coUapsing forward/ backward point if it is not L-collapsing for- 
ward/backward point for any L > 0; 

• oo-non-coUapsing point if it is a oo-non-collapsing forward and backward point 
simultaneously; 

• a collapsing point if it is either L-collapsing forward or backward for some L. 

Functions h'^ {ip , . The functions h~^{ip,ip') and h^{ip,%l}') are defined in M? as 
follows: h^{ip,ip') = Xq if the solution of Cauchy problem for Eq.Q with initial data 
ip{0) = ip, ipxif^) = "ip' collapses at value x = xq, xq > 0. By convention, we assume that 
h'^{'p,'p') = cxD if {'p,'p') is cxD-non-coUapsing forward point. Similarly, h~{'p,-p') = —xq 
if the solution of Cauchy problem for Eq.([6]) with initial data ^^{0) = -p, ipxiO) = 'p' 
collapses at value x = xq, xq < 0. 

The sets and Ul- We denote the set of all L-non-coUapsing forward points by 
and the set of all L- non-collapsing backward points by Uj^. In terms of the functions 

/i^ ('?/', '?/'') these sets are 

Ut = {(^,^0 e M^l h+iip,^') > L}, Ul = {{'Pj,,p') e M^l h-{ij,ij') > L}. 

The intersection of and U£ will be denoted by Ul- Evidently, if Li < L2 then 
W+ C Ul, Ul^ C and Ul, C Ul,. 

The values VL and Vl. We define 

= sup(u; — f/(x)), i7 = inf — f/(x)). 



2.2. Some statements about collapsing points 

In what follows we will use some statements from the paper [17], in particular so-called 
Comparison Lemma (reproduced in Appendix A for convenience). It is known [17] that 
for < Eq.([6]) has no bounded on R solutions, therefore we restrict our analysis by 
the case > 0. Also it is known ([17], Lemma 2) that all cx3-non-collapsing points 
for Eq.([6]) are situated in the strip —\/VL < 'p < VVt. Theorem [2ll below gives more 
detailed information about collapsing points for Eq.(j2]). 
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Theorem Let the potential U{x) he continuous and hounded on M. Then 

for each L there exist ipi and ip'j^ such that the set is situated in the rectangle 



The proof of Theorem [2ll is quite technical. We postponed it in Appendix A 
Another important statement is as follows: 

Theorem [2].2. Let the potential U{x) he continuous and hounded on M and 
h^{iljo,ipQ) = L < oo. Then h^lipjip') is a continuous function in some vicinity of 
the point {ipo, i'o)- 



The proof of Theorem [2l2 can be found in Appendix B It is worth commenting 
Theorem |2l2 as follows: 

(i) Analogous statement is valid for the function h~{ip^ip'). 

(ii) It follows from Theorem[2l2 that if the potential U{x) is continuous and bounded 
on M then Uf^ and Ul are open sets. The boundary of the set lA^ consists of continuous 
curves and corresponds to the level lines h~^{ip,il!') = L of the function h'^{ip,ip'). This 
boundary consists of the points ip') such that the solution of Eq.([6]) with initial data 
ipifS) = ip, ipx{0) = i^' satisfies one of the conditions 

lim iIj{x) = +00 or lim iIj{x) = —00. 

x—^L x—>L 

Correspondingly, the boundary of the set U£ is also continuous and consists of the points 
such that the solution of Eq.(l6]) with initial data ip{0) = ip, ipx{0) = '4'' satisfies 
similar conditions 

lim tp^x) = +00 or lim = —00. 

x—^—L x—^—L 

(iii) Theorem O 2 does not impose any restriction to the behavior of {ip , in 
a vicinity of a point {ipo,'^o) where h^^ipOyi^'o) = oo- practice, this behavior may be 
very complex, see Sect J2.3.2l 

The set of solutions for Eq. ([6]) that collapse at a given point x = xo can be described 
more precisely in terms of asymptotic expansions. 

Theorem [2l3. Let x = xq he an arhitrary fixed real. Assume that U{x) = oj — U{x) 
in a vicinity of x = Xq can he represented as follows 

U{x) = Uo + Ui6 + U26^ + + [6^] 

where 6 = x — Xq. Then the solutions of Eq.lj^ which satisfy the condition 

lim ijj{x) = +00 (7) 

x^xo 

ohey the asymptotic expansion 

Pi 

^{b) = ^ + + ^2(5^ + ^3(53 In \b\ + + ^4(5^ + o (5^) . (8) 


Here C G M a free parameter and 

6 ' 4 ' 5 ' 6 V 12 V 
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Figure 1. The point P, the manifold M+(xo), the curve a+(a;o)- 



Proof: The result follows from straightforward substitution of series ([H]) into Eq.(l6]). 

■ 

Theorem [213 should be commented as follows. 

(i) The free parameter C is "internal" parameter of continuous one-parameter set of 
solutions which tend to +oo at the point x = Xq. This situation can be illustrated by the 
following heuristic reasoning. Let P = {ipo, i/j'q) be a collapsing point and the solution 
ipi^x) of Cauchy problem for Eq.Q with initial data ipi^O) = ipo, ipxi^) = "^o collapses 
at X = Xq > 0. Then, generically ip{x) belongs to a continuous one-parameter set of 
solutions which also satisfy the condition and obey the expansion ([H]). In 3D space 
{x^ip^ip') this set generates 2D manifold M^{xq) (see Fig{T]). Intersection of M+(xo) 
with the plane (■?/', if)') at a; = includes the point P and is non-empty. Generically, this 
intersection in some vicinity of P is an ID curve which we denote a;+(xo). In the plane 
{ipjip') this curve corresponds to the level line h'^{ip,ilj') = xq. 

(ii) Since Eq.(l6]) is invariant with respect to the symmetry ip — )■ —tp, the solutions 
of Eq. ([6]) which satisfy the condition 

lim ipi^x) = —CO 

x—>-xo 

obey the same, up to sign, asymptotic expansion (j8]). The corresponding manifold 
M~(xo) and the curve a~{xQ) are defined in the same way. 
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Figure 2. The sets 14^ and 14^^ for Eq.® and the parameters w and A lying in the 
first (panels A and B) and the second (panels C and D) gaps, see SectlS] for detail. 
The parameters are: (A) w = 1, A -1; (B) uj ^ I, A ^ -3; (C) uj = 4, A = -4; (D) 
Lo = 4, A = —10. The sets were obtained numerically by scanning of initial data plane 
for Eq.®. The areas Ut: = n U~ are shown in dark. 



2.3. Example: the sets lAi^ andlAi for the cosine potential 

Let us now give now examples of the sets Uf^ for Eq.06]) in the case of the cosine potential 
(E]). Eq.(l6]) takes the form 

'^xx + {oj- A COS 2x)i) - = 0. (9) 

The sets Uf^ possess the following symmetry properties: 

1. Since the nonlinearity in Eq.(l9]) is odd, both the sets are symmetric in the plane 
{'il'yip') with respect to the origin; 

2. Since Eq.dH]) is invariant with respect to x- inversion the sets and are related 
to each other by the symmetry with respect to the axis ip and due to p.l, with 
respect to the axis ip' also. 




2.3.1. Ti -non-collapsing forward/backward points of Eq.^. The sets were found 
by thorough numerical scanning in the plane of initial data [tp, tp') (some details of 
numerical procedure can be found in SectJS]). The numerical study shows that for any 
values of parameters u and A the sets are infinite curvilinear strips. The typical 
shapes of the sets for Eq. ([9]) are shown in Figj2l The boundary of is represented 
by two continuous curves a^. The curve consists of such points {iPo,iP'q) that the 
solution ipl^x) of the Cauchy problem for Eq.Q with initial data ip{0) = ipo, ip'{0) = iJjq 
collapses at x = vr and \imx^.,^ip{x) = +00. At the curve a~ the solution iIj{x) of the 
corresponding Cauchy problem obeys the condition limx^T^ iIj{x) = —00. Similarly, the 
boundary of U~ is represented by two continuous curves (3^. The curves (3^ consist of 
points {ipQ, ip'o) such that the solution ipi^x) of the Cauchy problem for Eq.(l9]) with initial 
data ip{0) = ipo, ip'{0) = ip'^ collapses at x = — vr and \imx^_^ilj{x) = ±00. 

vr-non-collapsing forward anc? backward points of Eq.(l9]) form the set = W^flW^. 
It follows from Figj2] that this set may consist of several disjoined components. More 
detailed discussion of the sets and U-,, is postponed to SectlSl 
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2.3.2. TTn-non-collapsing forward/backward points of Eq.^, n> 1. FigJ3] exhibits the 
sets Uq^ for a; = 1 and various values of A. The sets W^r are the reflections of the sets 
14q^ with respect to the axis. It follows from Figj3] that the sets Uq^ have quite a 
complex layered structure. When n grows, the structure of U^^ becomes more complex 
resembling fractals. The situation is similar to one described in [19] for Eq.dl]) in the 
case of delta-comb potential. 

3. Symbolic dynamics: theory 

In this section we give a theoretical background for description of the non-collapsing 
solutions of Eq.(l6]) in terms of symbolic dynamics. Results of such kind are well-known 
in dynamical system theory. The language and the technique go back to 60-70-ties, 
see e.g. [20l EH |22]. In fact, the conditions which we formulate (Theorem |3ll) can be 
regarded as some version of the Conley-Moser conditions, see e.g. |22]. A peculiarity of 
the statement which we give below is that it is convenient for direct numerical check in 
practice. 

Let [ipjip') be Cartesian coordinates in and fi{S) be a measure of set S in M^. 
Remind that a function f{x) is called 'y-Lipschitz function if for any Xi and X2 the 
relation holds 

\fix2) - f{xi)\ < 71x2 - Xi\. 

Also introduce the following definitions. 

Definition. Let 7 be a fixed real. We call an island an open curvilinear quadrangle 
D G M."^ formed by nonintersecting curve segments a~^, /?"'", a~ , P~ {a^ and a~ are 
opposite sides of the quadrangle and have no common points as well as /3+ and f3~) such 
that 

• the segments a"*" and a~ are graphs of monotone non-decreasing/non- increasing 
7-Lipschitz functions ip' = v±{ip); 

• the segments and (3~ are graphs of monotone non-increasing/non-decreasing 
7-Lipschitz functions = w±{ip'); 

• if v±{tp) are non-decreasing functions then w±{tp') are non-increasing functions and 
vice versa. 

Definition. Let 7 be a fixed real and D be an island bounded by curve segments 
a"*", /3~^, a~ , /3~ . We call v-curve a segment of curve (3 with endpoints on a~ and a"*" 
which 

• is a graph of monotone non-decreasing/non-increasing 7-Lipschitz function ip' = 

• if are graphs of monotone non-decreasing functions, then v{iJj) is also a monotone 
non-decreasing function. If are graphs of monotone non- increasing functions, 
then v{ip) is also a monotone non- increasing function. 
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Figure 4. An island D with v-curve /3, v-strip V, h-curve a and li-strip H. 

Similarly, we call h-curve a segment of curve with endpoints on /3~ and which 

• is a graph of monotone non-increasing/non-decreasing 7-Lipschitz function = 

• if are graphs of monotone non-decreasing functions, then w{ip') is also a 
monotone non-decreasing function. If are graphs of monotone non- increasing 
functions, then w{i/j') is also a monotone non- increasing function. 

Definition. Let D be an island. We call v-strip a curvilinear strip contained between 
two nonintersecting v-curves, including both v-curves. Similarly, we call h-strip an 
open curvilinear strip contained between two nonintersecting h-curves, including both 
h-curves. 

FigJH illustrates schematically the definitions introduced above. 

Let us denote fl^ the set of bi-infinite sequences {. . . ,i_i,io,ii, . . .} where ik G 
{1,2,..., A^}. fi^ has the structure of topological space where the neighborhood of a 
point a* = {. . . ,i*_^, Iq, il, . . .} is defined by the sets 

Wk{a*) = {aeQ''\t,=t*,\j\<k}. 

Let T be a diffeomorphism defined on a set D = [J^^i Di where each Di C M^, 
i = 1, . . . ,N, is an island and all the islands are disjoined. Introduce the set V of 
bi-infinite sequences (called orbits) 

s = {. . . ,p_i,po,Pi, • • •}, Tpi = pi+i, 

where each pi = {ipi, ip'^), i = 0, ±1, ±2, . . ., belongs to D. V has the structure of metric 
space with norm defined as 
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Define a map S : P — as follows: ik is the number i of the component Di where 
the point pk lies. The following statement is valid: 

Theorem [3l 1 . Assume that: 

(i) a diffeomorphism T is defined on a set of N disjoined islands D^, i = 1, . . . , N , 
D = Uf=i A; 

(a) for any i, i = 1, . . . , N , and for each v-strip V E Di the intersection TV fl Dj, 
j = 1, . . . , N is non-empty and is also a v-strip. Similarly, for any i, i = 1, . . . , N , and 
for each h-strip H E Di the intersection T^^H fl Dj, j = 1, . . . , N is non-empty and is 
also an h-strip. 

(Hi) for the sequences of sets defined recurrently 

A,=D, A;^=T-'A-_,nD 

the conditions hold 

lim /i(A+) = 0, lim /i(A;) = 0. 

n— ^-oo n— >oo 

Then 11 is a homeomorphism between the topological spaces V and . 



Proof of Theorem Ol is postponed in Appendix C 



4. Coding of solutions 

4.1. Poincare map 

Assume now that the potential U{x) is continuous and ir-periodic 
U{x + n) = U{x). 

The Poincare map T : — )■ associated with Eq.(l6]) is defined as follows: if 
p = {ipyip') G then Tp = {ip{7[),'ipx{Tc)) where ip{x) is a solution of Eq.(|6]) with 
initial data ipi^O) = ip, 'V^x(O) = 'ip'. 

The map T is an area-preserving diffeomorphism. It is important that T is defined 
not in the whole M^, but only on the set of vr-non-collapsing forward points for Eq.([6]), 
i.e. W+. Inverse map T^^ is defined on the set W^. Evidently, for each p G W+ the 
image Tp G U~ and for each q G U~ the pre-image T^^q G W^j^, therefore = U~ 

and T-^U- = U+. 

If, in addition, the potential U{x) is even, U{x) = U{—x), Eq.(|6]) is reversible. 
The prototypical example is the cosine potential (j5]) which appears as a basic model in 
numerous studies. Denote / the reflection with respect to ip axis in the plane {1^,1^'). 
Due to reversibility of Eq.Q, if p G W+ then 

T-^Ip = ITp. (10) 

Therefore the sets and are connected by the relations JW^ = W^, JW^ = 
The set = nW~ consists of the points which have both T-image and T-pre-image. 



Coding of nonlinear states for the Gross- Pitaevskii equation 



12 



Theorem [21 1 implies that is bounded. It follows from Sect J2. 31 that may consist 
of several disjoined components Di, i = 1, . . . , N. 

The orbits defined by T are sequences of points (finite, infinite or bi-infinite) {pn}, 
such that Tpn = Pn+i- The fixed points of T correspond to vr-periodic solutions of 
Eq-© (such solutions do exist for quite general periodic potential U{x), see [23]). For 
a fixed point p let us denote DTp the operator of linearization of T at p. Let Ai^2 be the 
eigenvalues of DTp. Since the map T is area-preserving, A1A2 = 1. Depending on the 
behavior of T in a vicinity of a fixed point, it may be of elliptic or hyperbolic type |22j . 
In the case of hyperbolic fixed point both A1.2 are real and in the case of elliptic point 
they are complex conjugated, |Ai^2| = 1- Also we call a k-cycle an orbit which consists 
of points pi, . . . ,pk E M'^ such that 

Tpi=p2, Tp2=P3,---, Tpk=pi. 

Evidently pi,...,pk are fixed points for T'^. The A;-cycles correspond to /cvr-periodic 
solutions of Eq.([6]). A k-cjcle also may be of elliptic or hyperbolic type. This is 
determined by the type (elliptic or hj^erbolic) of the fixed point pi for the map T^. 

Below we consider bi-infinite orbits which lie completely within the set U^^. Basing 
on Theorem Ol we formulate necessary conditions which guarantee that these orbits 
can be coded unambiguously by the sequences of numbers i of in the order the orbit 
"visits" them. 

4-2. Symbolic dynamics: application to Eq.^) 

The application of Theorem [3ll to Eq.([6]) gives sufficient conditions for existence of 
coding homeomorphism. They can be formulated as follows: 

Hypothesis 1. The set consists of disjoined islands D^, i = 1, . . . , N, i.e. of 
curvilinear quadrangles bounded by curves which possess some monotonic properties 
(see the definitions in Sectj3]). 

Hypothesis 2. The Poincare map T associated with Eq.([6]) is such that 

(a) T maps v-strips of any Di, i = 1,. . . ,N, in such a way that for any v-strip V, 
V G Di, all the intersections TVCiDj, j = 1, . . . , N, are nonempty and are v-strips. 

(b) the inverse map maps h-strips of any Di, i = 1, . . . , N, in such a way that for 
any h-strip H, H E Di, the intersections T~^H fl Dj, j = 1, . . . , N, are nonempty 
and are h-strips. 

Hypothesis 3. The sequences of sets defined as follows 

A+=w., A+ = rA+_inw., 
A,=u^, A- = T-'A;^_,nK, 

are such that lim„_5.oo /^(^n ) = 0- 

It follows from Theorem Ol that if Hypotheses 1-3 hold then one can assert that 
there exists a homeomorphism between all bounded in M solutions of Eq.(l6]) and the 
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sequences from which can be regarded as codes for these solutions. The verification 
of Hypotheses 1-3 can be done numerically. From practical viewpoint, the following 
comments may be useful: 

1. If the periodic potential U{x) is even the point (b) of Hypothesis 2 follows from 
the point (a). In fact, if H is an h-strip then IH is a v-strip where / is a reflection with 
respect to ip axis. Then the statement (b) follows from the relation ( |T0|) . 

2. Let DTp be the operator of linearization of T at point p. Let 



(11) 





Define the functions 

gi{p) = (DTpei,ei) ■ (DTpe2,ei), 

g2{p) = pTpei,e2) ■ pTpe2,e2). 

Then the following statement is valid: 

Theorem 14]. 1. Assume that the potential U{x) is even and the following conditions 
hold: 

• is an infinite curvilinear strip; 

• W+ n U- =U^ = U r=i D^ where Di are non- overlapping islands; 

• for each pair {i,j), if 

— [5f are graphs of monotone non- decreasing functions then for any p G T^^DjH 
Di the relations gi{p) > 0, g2{p) > hold; 

— I3f are graphs of monotone non-increasing functions then for any p E T^^DjH 
Di the relations gi{p) < 0, g2{p) < hold. 

Then the conditions of Hypothesis 2 take place. 



The proof of Theorem HI 1 can be found in Appendix D It follows from Theorem 
mi that numerical evidence for Hypothesis 2 can be given by calculation of g\{p) and 
5^2 (p) within the set IA.„. 

3. If the periodic potential U{x) is even the relation ( JTOl) implies that yu(A+) = 
/i(A,^) for any n. Therefore in order to verify Hypothesis 3 in this case it is enough to 
check the condition lim„_>oo /^(-^^ ) — only. 

In the next section we describe the results of numerical study for Eq. ([6]) with cosine 
potential ([5]). We present numerical evidence that Hypotheses 1-3 hold for vast areas in 
the plane of parameters (a;. A). 



5. The case of cosine potential 



For numerical study of Eq.(l6]) with cosine potential ([5]), i.e. of Eq.Q, special interactive 
software was elaborated. It is aimed to fulfill thorough numerical scanning of the plane 
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{ip, i^') of initial data and visualize the sets Uj^ and Ul for a given L. Also the software 
allows to measure areas of Ul, to trace orbits generated by iterations of T, to find fixed 
points of T^, = 1, 2, . . ., to calculate values gi,2{p) ^^id visuahze areas where (71,2 (p) > 
and 5'i,2(p) < 0, and has some other useful features. 

For the numerical scanning of the plane {tp^ip') a 2D grid with steps Aip, Aip' was 
introduced. For each initial data -0(0) = ipQ, 'ip'if^) = ip'o at the grid nodes the Cauchy 
problem for Eq.Q on the interval [0; L] was solved numerically. If the solution tplx) 
of the Cauchy problem remains bounded (in modulus) by some large number B on the 
interval [0; L] we concluded that no collapse occurs and these initial data were regarded 
as an L-non-coUapsing point. Typically the values Atp = 0.0005 and Atp' = 0.0002 were 
taken. We found that the results for B = 100 and B = 1000 in all the cases were almost 
indistinguishable. 

Let us set forth the results of the numerical study of Eq.Q for each of the 
Hypotheses separately. 

5.1. Hypothesis 1 

Some examples of the sets W+ and IA~ are shown in Fig|2l In all the cases the sets 
and U~ are curvilinear strips. We found that this is a general feature of Eq.Q for all 
values of the parameters u and A that we considered. 

The shape of the strips may be quite complex and their intersection Ut, may 
consist of different number of disjoined sets. Since the strips lA^ and U~ are related 
to each other by symmetry with respect to the ip axis, the typical situation is that 
consists of several number of curvilinear deltoids (see Figl2], panels B and D) which are 
symmetrical with respect to or ip' axes. 

Figj5] shows the regions in the parameter plane (w, A) where such decomposition of 
takes place. Due to the symmetry 

A -A, X X + 7r/2 

the study has been restricted to the area A < 0. The zones corresponding to gaps 
and hands are also shown. Let us remind that the separation of the gap zones and the 
band zones in the parameter plane (w, A) is the key point for the theory of linearized 
(Mathieu) equation [2lj 

V'xx + {uj - A cos 2x)'4> = (12) 

If a point {u,A) belongs to a band, all the solutions of Eq. f|T2l) are bounded in R and 
if it belongs to a gap, all of them are unbounded. It is known that band/gap structure 
also plays an important role in the theory of nonlinear equation (Q (see e.g. [1]). In 
terms of the Poincare map T associated with Eq.Q, if the point (w. A) is situated in a 
band, then the origin 0(0, 0) is an elliptic fixed point for T, and if {u, A) lies in a gap 
then 0(0, 0) is a hyperbolic fixed point for T. 

In Figl5] two curves marked as = 3 and N = 5 are depicted. In the area above 
the curve N = 3 and below the curve N = 5 the set consists of three connected 



Coding of nonlinear states for the Gross- Pitaevskii equation 



15 



CO 



6 




- 2 



4 







6 



-12 



-8 



-4 







Figure 5. The plane of parameters (tu, A) with band and gap zones. The boundaries 
of the regions where iY^r is disjoined in three (the curve iV = 3) or five (the curve 
N — 5) connected components are shown. These components are islands (in the sense 
of Sectl3|) if uj and A belong to the dark zones between the marked curves and upper 
boundaries of the corresponding gaps. 



components, in the area above the curve = 5 it consists of five connected components 
etc. Due to Theorem [2l2 the boundary of each of the component is continuous, but 
a conclusion about monotonicity and Lipschitz properties of the boundaries should 
be made using numerical arguments. Our numerical study indicates that all these 
components are islands in the sense of Sectl3]in the areas (in dark) between the marked 
curves and upper boundaries of the gaps, see FigjSl 

We note that possible numbers of islands are related (indirectly) to numbers of fixed 
points of the Poincare map T. In its turn, the number of fixed points of T is determined 
by the number of band or gap where the point {u, A) is situated. More detailed analysis 
of these relations is an interesting issue for a further study. 

5.2. Hypothesis 2 

Since the potential U{x) is even and U:^ and U~ are infinite curvilinear strips, the 
verification of Hypothesis 2 can be fulfilled using Theorem Ull. To this end the 
calculation of the values of gi{p) and g2{p) was incorporated into the procedure of 
the numerical scanning 1 described above. To confirm the results we also used direct 
visualization of the vectors DTpei 2 for various p G U.,,. 

Let us describe in detail the case u = 1, A = —2 which is typical of the grey zone 
situated in the first gap, see Fig|5l The set consists of three islands Di, D2 and 
D3, see Figjni Their pre-images T~^Di, T^^D2 and T^^D^ intersect Di, D2 and D-^. 
Figl6] shows the signs of gi,2{p) for all the intersections T^^Di n Dj, i,j = 1,2,3. The 
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^ ^ 

Figure 6. The region of tlie plane {ip,tp'), —1.5 < i/j < 1.5, —1.5 < ijj < 1.5 with the 
sets and T~^U.^, lo — 1, A — —2. The set Ut^ consists of three islands £>i, D2 and 
D3, the boundaries [3^2 3 of the islands are marked by bold dash lines. The areas where 
gi{p) > and gi{p) < for p S T~^Di^2,3 are indicated in panel A. Corresponding 
areas for 52 (p) are shown in panel B. 

boundaries (3^2 3 of ^^e islands are marked by bold dash lines. In the islands Di and 
D3 the boundaries /Sf^ are graphs of increasing functions whereas for D2 the boundaries 
(3^ are graphs of decreasing functions. It follows from Fig|6] that the signs of 5'i,2(p) 
conform the Hypothesis 2. 

Overall, the numerical study shows that Hypothesis 2, as well as Hypothesis 1, 
holds for u and A lying in the dark areas in FigJSl 

5.3. Hypothesis 3 

The behavior of /x(A+), n = 1, 2, . . ., for various values of parameters uj and A was also 
studied numerically. Some of the results are depicted in FigJTl 

It follows from FigJT] that Hypothesis 3 is valid not for all the cases under 
consideration. A natural obstruction for Hypothesis 3 to hold is presence of elliptic 
fixed points or cycles. Due to KAM theory, in vicinity of an elliptic fixed point (or 
cycle) there exists a set of positive measure that consists of points which remain in this 
vicinity after any number of iterations of T. This means that Hypothesis 3 is not valid if 
the point (w. A) is situated in a band in the plane of parameters (see Sect JS.ip . because 
in this case the point 0(0,0) is an elliptic fixed point of T. This situation takes place 
for the case 1, w = 1 and A = —0.7, in FigJTl 

As the point (cu. A) crosses a lower boundary of a gap in the plane of parameters, 
the point 0(0, 0) becomes a hyperbohc fixed point and a pair of elliptic 2-cycles appears. 
It have been observed that this bifurcation is the first bifurcation in a cascade of period 
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10 20 30 n 10 20 30 

Figure 7. Panel (A): The area of the set A+ versus n for w = 1 and (1) A = —0.7, 

(2) A = -1.2 and (3) A = -2. Panel (B): The ratio i?„ = /i(A+_^J//i(A+) for the 
same values of u! and A. Only few points are shown in panel (B) for the cases (2) and 

(3) since accuracy of calculations drastically falls due to division of small numbers. 

doubling bifurcations. Each of the bifurcations of this cascade gives birth to elhptic 
cycles of double period. Omitting the details, we summarize that the gap zones in the 
plane {u, A) also contain areas where Hypothesis 3 is not valid due to presence of the 
elliptic cycles. 

At the same time, if the point {u, A) is situated in the dark zones in FigJSl numerical 
results show that Hypothesis 3 holds. Moreover, our results allow to suppose exponential 
convergence of /i(A+) to zero. The ratios i?„ = /^(An+i)//^(^n) shown in panel B 
of FigJTl For the cases 2 and 3, these ratios are smaller then 1 and remain close to the 
value /i (TUn nUj,)/ fi (U-,, ) . 

To summarize, basing on the numerical results presented above one can conclude 
that if oj and A are situated in the dark zones of the parameter plane, see Fig\^ the 
conditions of Hypotheses 1-3 hold. Therefore for these values of parameters all the 
nonlinear states of GPE with cosine potential can be put in one-to-one correspondence 
with codes from Vt^ . When crossing the lower boundary of the grey zones (marked 
iV = 3 or iV = 5 in Figj5]) the conditions of Hypotheses 1 fail whereas other two 
Hypotheses remain valid. 

5.4- More details about the map T 

Let us describe in more detail the transformation of the sets and by the map T. 
The action of T on vertical sections of is shown in FigJHl Note that TU^ = 
= /W+ where / is the reflection with respect to ip axis. Fig J9] illustrates the mapping 
oiU^ and the islands Di, D2, -D3 (the shape of U:^ was calculated for u = 1, A = —3). 
It is practical to represent the map T as a composition of three transformations: 




U:^ — 7- — 7- 52 — 7- W^, where 5*1 is an infinite horizontal strip and 5*2 is an infinite 
vertical strip, see FiglHl The "boundaries" of both, and Si, are in infinity. The 
transformation — )■ 5*1 is a deformation. The transformation of Si to 5*2 consists in 
stretching of 5*1 in one dimension and contraction in another in such a way that the 
boundaries transform into vertical lines but go to infinity. The transformation 
5*2 —7- is again a deformation. As a result, T maps the islands Di, D2, -D3 into 
infinite curvilinear strips. Each of these strips crosses the islands Di, D2, -D3 and each 
of the intersections is a v-strip. 

An interesting issue is the study of ordering of the v-strips in and the h-strips 
in corresponding to codes with coinciding blocks. The understanding of the strip 
ordering is also of practical usage, since it explains the order of the nonlinear modes as 
they appear in shooting procedure, see Sect J5.5[ We describe the ordering of v-strips, 
the ordering of h-strips is similar. Assume that U^^ consists of disjoined islands and 
A^ is odd. Consider orbits that visit the islands -Di_„, -Di_„+i5 • • • , i^i the given order. 
The points in Di^ which has this "prehistory" are situated in a strip Vi_^i_^^^-^^,,,i_2i-iio 
constructed by the following recurrence rule 

Vi i ^ = TDi n Di 

^_„i_„+i..i_2j-iio ~ ^^-ni-n+i---i-2«-i l~'-Dio- 



and (see [Appendix CD 



. . . C Vi_^j_^^j...i_2i_ijo C Vi_^_^^...i_2i_iio C . . . C Vi_2i_-^ig C Vi_ji(, C -Dip 
The orbit of a point p G V^_„i„„+i..i_2i_ijo has in its code a block 
(. . .i_„i-„+i . . .i-2i-i«o • • •)• 
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Figure 9. The action of tlie map T on and tlie islands Di, D2, D3 {the islands 
and their images are indicated by the same letters). The map T can be regarded as 
a composition of three transformations: W+ Si ^ S2 ^ U~: (i) deformation of 
W+ into infinite strip 5*1, the "boundaries" (5"^ of both, and 6*1, are in infinity; (ii) 
stretching of Si in one dimension and contraction in another one to get an infinite strip 
52, the boundaries become vertical lines but go to infinity; (iii) deformation of 
^2 into U~ . As a result, T maps the islands D2, -D3 into infinite curvilinear strips. 
Each of these strips crosses again the islands D2, D3 and each of the intersections 
is a v-strip. The shape of W+ was calculated for = 1, A = —3. 

Since each of ik can take the values 1, . . . ,N, there are A^^^^ strips in each coded by 
the sequence of length n + l. The algorithm for their ordering in follows immediately 
from geometrical properties of the intersection of the strips U~ and U^. It can be 
described as follows: 

1. Mark the islands Di, . . ., Dj^ as they are ordered in see FigHOl Draw an 
arrow Jq over all of them pointing from to Dat. 

2. Draw arrows Iqi, . . . , Iqn over each island in such a way that the directions 
of rightmost Jqi and leftmost Jqat arrows coincided with the direction of /q, but the 
directions of any two neighboring arrows were opposite. Sketch v-strips in Di^ in such a 
way that their ordering (from Vu^ to V^ig) agreed with the direction of the arrow Joig. 

3. Draw an arrows over each of ^.^jg by the same manner and sketch v-strips 
^_2i_iio according the directions of these arrows, etc. 
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Figure 10. Algorithm for ordering of v-strips in . 



It turns out that the algorithm given above is similar to procedure of ordering of 
localized modes for DNLS described in \25\ . 

5.5. Examples 

Let the parameters (cu, A) belong to the dark zone in the first gap, see FiglS] Then 
all the bounded in R solutions of Eq.(|9]) can be coded by bi-infinite sequences of three 
symbol "alphabet". Conversely, for each bi-infinite "word" composed by symbols of 
this "alphabet" these exists a solution with corresponding code. The symbols may be 
chosen as "— ", "0" and "+", and they mark entering of the orbit into Di, D2 and 
respectively. 

Example 1. Localized modes described by Eq.(|9]) correspond to the codes with 
finite numbers of nonzero symbols. In particular, Eq.(|9]) admits well-known solution in 
the form of bright gap-soliton, ip{x), [71 19], locahzed in one well of the potential, see 
Fig{TT] A. This solution corresponds to the code (. . . 00 + 00 . . .). Also there exist the 
gap soliton solution —ipi^x) with the code (. . . 00 — 00 . . .). 

Example 2. There exist exactly two vr-periodic solutions of Eq.(l9]) with the codes 

(... + + + ...) and (. . . . . .), related to each other by symmetry ip — )■ —tp, see 

FigilllB. 

Example 3. There exists dark soliton solution of Eq.(|9]) corresponding to the code 

(. . . h + + . . .), see FiglH] C. Also the coding predicts that there exist other 

solutions of this type, having the codes (. .. + + + ...), (. .. 00 + + + ...), 

etc. 

Example 4- There exist "domain wall" -type solutions corresponding to the codes 

(. . . 000 + + + ...), (. . . 000 . . .). These objects were found to exist in the case of 

GPE with attractive nonlinearity |12]. The coding approach predicts their existence in 
the case of repulsive nonlinearity also. They have been found numerically, see Fig. [H] 
D. 

Example 5. Consider boundary value problem for Eq.Q on the interval [— 47r, 0] 
with Neumann boundary conditions at x = —An and x = 0. These solutions can be 
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Figure 11. Nonlinear states for GPE, u = 1, A = —2: A: briglit gap-soliton, tlie code 
(. . . 00 + 00 . . .); B: 7r-periodic structure, tlie code (... + + + ...); C: dark soliton, the 
code (. . . h + + ...); D: domain wall, the code (. . . 000 + + + . . .)• 

viewed as reductions to the interval of length An of periodic solutions with period Stt 
which satisfy additional symmetry conditions 

u{x) = u{—x); u(— 477 + x) = u(— 47r — x). 

The codes for these solutions are of the form 

(. . . 6*56*4 6'3^2 ^1^2^3^4 6*56*4 6*3^2 , ^1^2^3^4^5 • • •)) 

the period 

where 6i, i = 1-^5, is one of the symbols "+", "0" or "— ". Therefore there are 3^ = 243 
solutions of this type. In FiglT2]nine of these solutions are depicted for A = —2, u = 1. 
From numerical viewpoint, these nonlinear modes can be found by shooting method 
taking initial data i(j{—4:7r) = ip, i/jx{—4:7r) = and adjusting ip in such a way that 
■^xlO) = 0. As ip increases, one gets the nonlinear modes one by one in the order 
described in Sect J5.4[ 

All the solutions shown in FiglT2] have the codes with ^1 = 6^5 = "0" . Therefore 
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they can be viewed as approximations for localized modes which have the domain of 
localization of length Att. These localized modes correspond to the codes 

(...00 ^2^3^4 00...). 

There are 3^ = 27 sequences of this type but only 10 of them (including one which 
consists of zeroes only and corresponds to the zero solution) are different in the sense 
that they are not related to each other by symmetry reductions. FiglT2] shows just these 
nine nonzero solutions. 

6. Conclusion 

In this paper, we describe the method for coding of nonlinear states covered by ID 
Gross-Pitaevskii equation with periodic potential U{x) and repulsive nonlinearity. We 
prove that under certain conditions there exists one-to-one correspondence between the 
set of all bounded in R solutions of Eq. ([6]) and the set of bi-infinite sequences of numbers 
1, . . . , A^. These sequences can be regarded as codes for the solutions of Eq.(j6]). The 
number is determined by the parameters of Eq.(l6]). It is important that (i) each 
coding sequence corresponds to one and only one solution and (ii) each solutions has 
corresponding code. The conditions for the coding to be possible are presented in a form 
of three hypotheses. For a given U{x), the hypotheses should be verified numerically. 
We report on numerical check of these hypotheses for the case of cosine potential, i.e., 
for Eq.Q, and indicate the regions in the plane of parameters {u, A) where the coding 
is possible. 

Heuristically, the coding technique described above can be interpreted as follows. 
The periodic potential can be regarded as an infinite chain of equidistantly spaced 
potential wells. It is known that if U (x) is a "deep enough" single- well potential, Eq. ([6]) 
admits one or more localized solutions called "fundamental gap soliton", FGS, in [THITS] . 
Also Eq.(|6]) admits zero state. Assume that in total there exist N states (including 
zero state) described by the single well potential. Note that the number is odd: 
since the nonlinearity is odd if ipi^) is a FGS then —ip{x) is also a solution Eq.([6]) 
and the zero solution should also be taken into account. In these terms the coding 
means that one assigns to each of possible single-well states a number from 1 to and 
attributes to "bound states" of these entities situated in the wells of periodic potential 
bi-infinite "words" consisting of numbers from 1 to A^. This viewpoint exploits an 
analogy between periodic problem and discrete problem replacing the solution on each 
period by a lattice node. The corresponding reduction can be made consistently and 
rigorously using Wannier functions technique [26] but the resulting system of discrete 
equations is nonlocal and quite difficult for a comprehensive study. 

The approach presented in this paper may be applied to Eq. ([6]) with different types 
of the periodic potential U{x). Also it may be extended in various directions. In 
particular, preliminary studies show that it can be applied with minor modifications to 
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Figure 12. The solutions of Eq.® = -2, w = 1) on the interval [-47r, 0] with 
Neumann boundary conditions at x = —A-k and x — Q (nine panels below) and the 
position of corresponding points in the intersection of U^^ and U^^ (upper panel). The 
upper panel represents the scaled rectangle marked in panel D of Fig|3]with fragments 
of U^^ (in black) and U^^ (in grey). The codes (the blocks {O162636465)) are: A: 
(0 + -00); B: (0 + - + 0); C: (0 + 000); D: (0 + + 0); E: (0 + - 0); F: (0 -f -f 00); 
G: (0 + + 0); H: (0 - 0); I: (00 + 00). 
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the equation 

i^cccc + {uj-U (x))^ + g,ij' + g^ij' = 0. (13) 

where g2 < 0. Eg. (1131) also arises in the theory of BEC [271 128] . The shapes of the sets 
in this case are similar to ones described above for Eq. ([6]) . Another possible extension of 
this approach can be made for comj^leo; nonlinear states of GPE, ip{x) = p{x) exp{i(j){x)}. 
It is known [9j that the amplitude p{x) obeys the equation 

Pcccc + {oJ- U{x))p - p3 _ = 0, (14) 

p^ 

where M is an arbitrary real constant. For a given amplitude p(x), the phase (/)(x) can 
be found from the relation p'^{x)(j)x{x) = M. Other possible extensions of the approach 
may be related to the cases when U^^ consists of partially overlapping islands or of more 
general sets which are not islands at all. 

Having in hand a complete description of nonlinear modes for Eq.(l6]) in terms of 
their codes, one can return to the problem of stability of these modes. In our opinion, 
a relation between the code and the stability of corresponding mode is an interesting 
issue for further study. A good example of such a study is paper |29j where the similar 
problem was considered for DNLS. 

At last, let us note that the approach developed in this paper cannot be applied (at 
least, directly) to the case of GPE with attractive interactions, a < in Eq.(jl]). One 
can prove that in this case all the solutions of Eq.(jl]) are non-collapsing under quite 
general assumptions for the potential U{x). 
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Appendix A. Proof of Theorem Ol 

The following statement proved in [17J will be used below: 

Comparison Lemma. Let the functions y{t) and x{t), t G [a; 6] he solutions of 
equations 

yu-g{t,y) = (A.l) 

xtt-fit,x) = (A.2) 

correspondingly. Let also the following conditions hold: 

(i) f{t,^), g{t,^) are defined on [a; b] x [A; B] and locally Lipschitz continuous with 
respect to ^, ^ G [A; B], (A,B, b may be finite or infinite); 
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(^V 9it, > fit, for any te[a;b],^e [A; B]; 

(Hi) /(t, monotone nondecreasing with respect to ^, ^ G 

Let A < x{a) < y{a) < B and Xt{a) < yt{a). Then Xtit) < ytif) and x(t) < y(t) 
while A < x(t),y{t) < B or for the whole interval t G [a; h] 

In what follows we assume that the potential U (x) be continuous and bounded on 
M and use the notations introduced in Sect J2.1[ To prove Theorem Ol we need the 
following lemmas: 

Lemma [2ll. Let -i/^o > 0. Then for ipo > 

h^i^o, i^'o) < HiiJo) ^ r I "^t! (A-3) 

and Hq^iPo) ~ ^''X^"' ''J^hen ipo — )■ +oo. Here K(-) is complete elliptic integral of the 
first kind. 

Proof: Consider the equation 

0,, + n0 - 03 = (A.4) 

The solution (j){x) for Eg. flA.4p with initial data 0(0) = i/jq, 0^(0) = ■00 can be written 
in implicit form as follows 

X = ^=^fi==; A = 2(^;)2 + 2n^o'-^o- (A.5) 



The solution 0(x) tends to +oo at the point 

V2 dr] 



Xo 



and the following estimation holds 

, N r V2dr] K(l/V2) 

J^O ^^4 _ ^4 _ 2^(^2 _ ^2) ^0 

Now let us consider the solution ip{x) of Eq.([n]) with initial data -0(0) = -00; i'xiO) = "00- 
Since for ^ > y/Ti the function F(^) = ^^ — il^ is monotonic and ^^ — {cu — U (x))^ > F{^) 
one can apply Comparison Lemma from [17J to Eq.([n]) and Eq.f lA.4p . Therefore for 
X > the inequality 0(x) > 0(x) holds. This means that ip{x) collapses at a point 
/i'*'(0O)0o) ^ Xq < Hq^iPq). This proves Lemma [2ll. ■ 

Lemma [2l2. For each L there exists a value ipL such that the setUi is situated in 
the plane = (0,0') in the strip —ipL < < '0l- 

Proof: Due to Lemma Ol for each L there exists ipi such that there are no points of 
Ul in the sector > ip^, ip' > 0. Since Eq.([6]) is invariant with respect to the symmetry 
— )• —-0 the estimation (IA.3|) holds also for 0q < and 0o < — V^, therefore there are 
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no points of Ul in the sector ip < —ipL, V^' < 0. Making the transformation x — )■ —x 
and repeating the reasoning of Lemma [21 1 we obtain the estimation 

for the two cases: (i) ipQ > V^l, ipQ < and (ii) ipo < — V^, ipQ > 0. Similarly, 
hoii'o) ~ ~^ — oo. Therefore there are no points of Ul in the sectors 

ip > ipL, V^' < and ip < —ipL, 4'' ^ 0- This implies the statement of Lemma [2l2. ■ 

Proof of Theorem\^l: Due to Lemma [21 2 there exists the value ipi such that no 
points oi Ul are situated out of the strip —ipL < i' < ipi- Therefore it is enough to 
prove that there are no points of in two half-strips 

S^, = {(^,^0 G -^z. < ^ < ^L, ^' > R} 

= G M^l - < < ^L, ^' < -R} 

for R large enough. Let us prove this fact for S"^, the proof for S]^ is analogous. It 
follows from Lemma [2l2 that there exists the value ipL/2 > V'l such that all the points 
('0,V'') for > 4'L/2 and ip' > are L/2-collapsing forward points. Introduce the value 

Ml = min ^ (^^ - r/^) 

Evidently Ml < 0. Consider the solution tp{x) of Cauchy problem for Eq.(l6]) with initial 
data V(0) = V'o, V'o e [-^l] ^Pl] and tp^{0) = 'ipQ, ^'q > R. While -i'L/2 < ^{x) < i'L/2 
one has ipxxi^) ^ Ml- Then the following relations hold 

Mx) > ^x(O) + Mlx >R + Mlx, (A.6) 

ij{x) > ipiO) + Rx + > -ipL + Rx+ ^x^- (A.7) 

We claim that if the initial data for Eq.(|6|) are situated in the half-strip with 

R>ij'L = max |-^, ^ (8^L/2 + - M^L^) | 

then the segment of curve {{ip{x),ilJx{x)), < x < L/2} crosses the line ip = i^L/2 in 
some point where ip' > 0. In fact, assuming that —ipL/2 < ip{x) < ipL/2 for x G [0; L/2] 
from ( 1A.6I) and ( 1A.7I) one concludes that 

ipx{x) >R + MlX > 0, a; G [0; L/2]; 

i.e. we arrive at the contradiction. Therefore there exists a value x G [0; L/2\ such that 
ip{x) > 4'l/2 and '^/'^^(x) > i.e. {■ilj{x),ipx{x)) is L/2-collapsing forward point. Then 
{ip{0),tpx{0)) is L-collapsing forward point. So, for R> ip'^ there are no points from Ul 
in 5+. ■ 
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Appendix B. Proof of Theorem [2].2 

Proof of Theorem{^2. Introduce the following functions: 

(a) the function H~^{ip,'ip' ,Q;t) defined as follows: H~^{'ip,'ip' ,fl;t) = xq if the 
solution of Cauchy problem for the equation 

(j)^^ + fi0 - 0^ = 

with initial data 0(t) = ip, (f)x{t) = ip' collapses at the value a; = Xq, Xq > 0. Exact 
formula for H~^{tp, ip', Q; t) is 

H-^i^P, n;t)=t + r ^"^"^ (B.l) 

^2 



It follows from fIB.ip that if H~^{ipo,tpQ,QQ]t) < oo then for t fixed the function 
H^{ip,ip' ,Q;t) is a continuous function of the variables ip,ip',Q in some vicinity of 
the point (ipoji'o^^o)- 

(b) the function {ip , ip' ; t) defined as follows: h'^ {ip , ip' ; t) = Xq if the solution of 
Cauchy problem for Eq.([6]) with initial data ipit) = ip, ipxit) = %p' collapses at value 
X = xq, xq > 0. Evidently, if ip{x) is a solution of Eq.(l6]) then 

h+{^{o),Mo)) = h^{m,Mty,t) +t (B.2) 

(c) the two functions 

Qi(t) = min (cj — f/(x)), Q2{t) = max (w — ?7(x)) 

xe[t;L] xe[t;L] 

which are continuous functions in some vicinity of the point t = L. 

Also let us denote by Ds{C, (') a disc in with center at {(, (') and radius 6. 
It follows from the conditions of Theorem [21 2 that the solution ip{x) of Eq.([n]) with 
initial data ip{0) = ipQ, V'x(O) = ipQ satisfies one of the conditions 

lim tp{x) = +00, or lim ip{x) = — oo 

x-^L x—^L 

Let the behavior of ip{x) in vicinity of x = L obey the first of the two formulas above 
(the analysis of the second case is similar). Then there exists t such that ip{x) > ^/Tl^Jt) 



and ipx{x) > for x E [t;L). By virtue of Comparison Lemma, see Appendix A 
has 



one 



H^i^,^',n^{t*);t*) < h+i^,ij';t*) < fi2(r); t*) (B.3) 

for any t* G [t; L) and for any tp,ip' in some vicinity of the point {ip{t*),ip'{t*)). 

Let us describe 3-step algorithm which allows by given e > to find 6 > such 
that if {ip,ip') e Dsiipo,^) then 

\h+{ij,ij')-h+{iJo,ij'o)\<e. (B.4) 

1. By given e one can find t* such that the inequality holds 

^'(t*), Q2{t*); t*) - H+m*),^'{t*), Q,{t*); t*)\ < e/2 (B.5) 
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2. Since if+ (■?/', t/^', fi; t) is continuous there exists 6i > such that when {ipjip') G 
Dsi{'ipit*),ip'(t*)) the inequahties hold 

\H+i^, fi2(r); t*) - H^m*), ^'(r), fi2(r); t*)\ < e/2 

|#+(^, fii(r); f) - H+m*),^'{t*), n,it*y, t*)\ < 6/2 

It follows from (El]) and (EJ) that if (^,^') G ^^^(^(t*), ^'(t*)) then 

(b.6) 

3. The flow defined by Eq.(l6]) generates a diffeomorphism Tt* : — )■ which 
maps a point {ip{t*),ipxit*)) to a point {■ilj{0),ipx{0)) where ipi^) is a solution of Eq.(l6]). 
Then there exists 6 such that Dsiipo^^'o) C T^.D^, ^^(r)). By means of (IRJ) 
and (1B.6P one concludes that for {ip,ip') G Ds{iPo,iP'q) the relation (1B.4I1 holds. Theorem 
[2j2 is proved. ■ 



Appendix C. Proof of Theorem [3ll. 

Before proving Theorem [31 1 we prove the following lemma: 
Lemma [3ll. Let D be an island. 

(i) Let D D D V2 D . . . he an infinite sequence of nested v-strips such that 

lim ^I{Vn) = 0. (C.l) 

n—^oo 

Then the intersection 

00 

n=l 

is a v-curve. 

(a) Let D D Hi D H2 D . . . is an infinite sequence of nested h-strips and 
lim fi{Hn) = 0. 

n— >oo 

Then the intersection 

00 

n=l 

is an h-curve. 

Proof of Lemma 1^1. Let us prove the point (i), the point (ii) can be proved 
similarly. Denote the v-curves which bound the strip Vn by f3^ (which lies closer to 
and (which lies closer to f3~^). Let the endpoints of be (situated at a") and 
B:^ (situated at a"*"). Let the endpoints of /3~ be A~ (situated at a~) and B~ (situated 
at a"), see Fig JCll 

First, we show that A~ — )■ A^ and — )• B^ as n — > 00. The sequence of 
points Ai,A2, ... is situated on the curve a~ from one side from the point A~l and 
is "monotonic" in the sense that for any n the point A~_^^ is situated on a" between 
the points A~ and A'^. Therefore it has a limit point A^ G a~ . The sequence of 




Figure CI. Illustration to the proof of Lemma |3l 1 . 



points Af,A2, ... is situated on the curve a~ from one side of the point and has 
similar monotonic property. Therefore it also has a limit point A^ G a~ . Suppose 
that y4+ 7^ A^. Then, since are graphs of monotone non- increasing/non-decreasing 
7-Lipschitz functions and a~ is a graph of monotone non-decreasing/non-increasing 
7-Lipschitz function, the area of Vn cannot tend to zero as n oo. Therefore 
y4+ = A^ = Aoo- In the same way one can introduce the limit points G a"*" 
and G and conclude that B^ = B^ = B^o- 

Second, let coordinates of A^o be {%Ija,'>P'a) ^'^d coordinates of B^o be {iPb,4'b)- 
Consider a real value z situated between ipA and ipB- Since A~ — )■ A^ and B~ — )■ 5+ as 
n — )■ cxD there exists rio such that for n > tiq both /3+ and (3~ intersect the line ip = z. 
Denote the points of intersections of and /3~ with the line ip = z correspondingly 
C^{z) and C~{z). Evidently, both the sequences {C^{z)} and {C^{z)} have limits as 
n —7- oo. Denote these limits C^{z) and C^{z) correspondingly. 

Assume that at some z = z* one has \C^{z) — C~{z)\ = Sc > 0. The relation 
C~^{z) 7^ C~{z) cannot hold in some vicinity of the point z = z*, otherwise there exists 
a set of nonzero measure which belongs to all the nested strips Vn and fiiVn) does not 
tend to zero as n — )■ oo. Therefore for any e there exist a value Zi, such that |zi — z*| < e 
and C~^{zi) = C~{zi). This means that (al least) one of the ratios 

\C+{z,)-C+iz*)\ ^ \C-iz,)-C-iz*)\ 
\zi — z*\ ' — z*| 

is greater than 6c /2e. Since e can be taken arbitrarily small, this contradicts the 
condition that are graphs of monotone non- increasing/non-decreasing 7-Lipschitz 
functions. This implies that C~^{z) = C~{z) = C{z) for all z, tpA < z < ipB- 

Third, each of the curves f5~ is a graph of a monotone non-increasing/non-decreasing 
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7-Lipschitz function. Passing to the limit n — oo we obtain a curve consisting of the 
points C{z), ipA < z < ipB- This curve is also a graph of monotone non-increasing/non- 
decreasing 7-Lipschitz function with the same 7, (see [22], Sect 4.3), i.e. v-curve.B 

Proof of TheoremlM 1 ■ Evidently, for each p G P the image s = Ep G fl'^ is defined 
uniquely. Let us prove that for each s G there exists unique p G P such that s = Ep. 
Consider a sequence s = {. . . , i-i,io, ii, . . .}, ik G {1, . . . , A^}. Let us find the location 
of the points p G Di^ such that T~^p G -Di_i, T~^p G Di_^ etc. It is easy to check that: 

- the points p G Di^ such that T~^p G Di_^ are situated in the set V^_iio = 
TDi_^ H-Di;,. Due to the condition (i) of theorem, Vi_-^jQ is a v-strip. Moreover, Vj^q and 
VkiQ have no common points if j 7^ k. 

- the points q G Di such that T~^q G -Di_2 are situated in the set V^_2i_i = 
TDi 2 n which is a v-strip. The points p G -Djo such that T~^p = g G -Dj_i, 
T~^p = T'-'^g G -Di_2 are situated in the set Vi_2i_jjQ = Tyj ^j-i ^ Ao which is also 
v-strip. Evidently, Vi_2i_iio C Vi_^iQ. Also, Vji_^^^ and Vki_-^ig have no common points if 

Continuing the process, we have nested sequence of v-strips 

• • • ^_(„+i)i_„...i_2j-ijo ^_„..i_2j-iio C . . . C Vi_2i_ijo C Vi_-^i„ C -Dig 

such that Vi_j^^jjj„,j...j_2j_ijo ~ ^^_(„+i)i-n---j-2i-i -^«o- Since 

/i(^i_„...i_2i-iio) < /^(^n) 

the area of v-strip Vi_„..i_2i_iio tends to zero as ri — )■ 00. According to Lemma Ol the 
intersection of these nested strips, Voo, exists and is a v-curve. 

In the same manner the nested sequence of h-strips can be constructed, 

where Hi^i^^^^i^i^^^ = T''^Hi^i^,,,i^ fl Ao- The area of the strip Hi^i^^^^i^ tends to zero as 
n — )■ 00 so according to Lemma Ol the intersection of these nested strips, ifoo, exists 
and is an h- curve. 

The orbit p E V corresponding to bi-infinite sequence s = {. . . ,z_i,'io,^i, • • •} is 
generated by T- and T~ ^-iterations of the intersection Hoo H V^o which according to the 
definition of h- and v-curves consists of one point. Therefore p exists and is unique. 

The continuity of S and follows from the following observations: 

- since T is continuous, if p*^^^ = {. . . ,p^}l,pl^\p^^\ . . .} and p*^^^ = 
{. . . ,p^^l,p^Q\pf\ . . .} are close enough in V (i.e. the points p^^'' and p^q'' are close 
in R^), then their S-images share the same central block \i\ < k for some k. Therefore 
they are also close in fi^-topology; 

- if s*^^) = Sp(^) and s*^^^ = Sp*^^) share the same central block \i\ < k for some k, 
the points Pq^'' and p'^^ are situated in the curvilinear quadrangle Vi_j....j_2i_ijo (iHigi-^ i^, 
so p*^^-* and p*^^^ are close in P-topology. ■ 
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= {a| 


a 


= xei + ye2, 


X 


>0,y>0}, 




= {a| 


a 


= xei + ye2, 


X 


>0,y>0}, 




= {a| 


a 


= xei - 


\-ye2, 


X 


>0,y<0}, 




= {a| 


a 


= XBi - 




X 


>0,y<0}, 




= {a| 


a 


= xei - 


l-ye2, 


X 


>0,y<0}, 


wl_ 


= {a| 


a 


= XBi - 


l-ye2, 


X 


>0,y<0}, 




= {a| 


a 


= XBi - 


l-ye2, 


X 


<0,y<0}, 




= {a| 


a 


= xei - 


^Z/e2, 


X 


<0,y<0}. 



Lemma [41,1. Let S G M."^ be a compact connected set, T be a diffeomorphism 
defined on S, the operator DTp be nondegenerate for all p E S and ei,e2 be the unit 
vectors defined by IHl]) . 

I. Let for all p G S the relations gi{p) > and g2{p) > hold. Then 

a. For all p G S one and only one of the following alternative conditions holds 



[Al] DT,[rI^] C 
[A3] DT^IrI^] C 



[A2] DT, 
[A4] DT, 



r™2 



C 



C 



b. there exists 7 > such that for any two points pi = {ipijip'-^^) G S and 



P2 

Tp2 



(^2,^2) G S , tpi < il)2 and ip'i < tp2 the images Tpi 



Q2 



•^2,(P2, 



are such that 



in the case [Al] 
in the case [A2] 
in the case [A3] 
in the case [A4] 



< 02 - 01 < 7(02 - 0i) 
< 01 - 02 < 7(01 - 02) 

< 01 - 02 < 7(02 - 0'l) 
< 02 - 01 < 7(0'l - 02) 



//. Let for all p G S the relations gi{p) < and g2{p) < hold. Then 

a. For all p G S one and only one of the following alternative conditions holds 



[Bl] DT^ 



r™2 



c 



c 



r™2 



H+1 

d2 . 



[B2] DTp[R_^\ C 
[B4] DTp\R^_J C 



[B3] DTp[R_^, ^ ^.pL 

b. There exists 7 > such that for any two points pi = [ipijip'i) G S and 
P2 = ("02, "02) G S, ipi < -02 and i[)'i > -02 the images Tpi = qi = (0i,0']^), 
Tp2 = q2 = (02, 02) are such that 



in the case [Bl] 
in the case [B2] 
in the case [B3] 
in the case [B4] 



< 02 - 01 < 7(02 - 0i) 
< 01 - 02 < 7(0'l - 02) 
< 01 - 02 < 7(02 - 0'l) 
< 02 - 01 < 7(01 - 02) 
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Proof of Lemma^ 1 . Let us prove the point la, the point Ila can be proved similarly. 
Evidently, the relations gi{p) > and (72 (p) > mean that both the vectors DTpBi, 
DTpe2 are situated in the same quadrant, or for any p E S. Let for 

two points pi and p2 these quadrants are different. Connect these points by continuous 
curve ( G S. Since DTpBi, DTpe2 depend continuously on p there exist a point p* G (, 
such that one of {DTpBi, e^), i = 1, 2, j = 1, 2 vanishes therefore gi{p*) = or g2{p*) = 0. 
This implies that the point I is valid. 

Let us prove the point lb. Assume that gi{p) > and (72 (p) > hold and the 
situation [Al] takes place, the situations [A2]-[A4] can be treated similarly. It follows 
from the condition [Al] and compactness of 5* that there exists a supremum 

7 = sup 1^, I 1 = DTpZ, z G M^^, p e S. 




Let ip2 > i'l and ip2 > i)[ and pi = P2 = i'4'2,i'2), Qi = Tpi = 

q2 = Tp2 = (02, 02)- Then 



1p2 - V^l 

where ||r||/||p2 ^ when ||j92 ~^ 0. This implies that for close enough pi 

and P2 one can choose 7 > 7 such that for corresponding qi and q2 

< - 0'l < 7(02 - 0l) (D.l) 

This ordering is transitive: from the relation (ID.ip and the relation < 03 — 02 < 
7(03 — 02) it follows that < 03 — < 7(03 — 0i). Therefore we can omit the words 
"for close enough" above and state that fID.ip holds for any pi^2 such that -02 > ipi and 
1^2 > ipi- So, the point lb under the assumption [Al] is proved. In the same manner 
the point lb can be proved for other three cases, [A2]-[A4]. The proof of the point lib 
consists in considering in the same manner the situations [B1]-[B4]. ■ 

Proof of Theorem^ 1. Since Eq.([6]) is invariant with respect to x- inversion, the 
strip is symmetric with respect to the origin. If a point (-0; ip') is situated on one 
edge, a^, of the strip lA^ the solution ip{x) of Cauchy problem with initial data (-0; ip') 
obeys the condition 

lim'0(a;) = +00 (D.2) 

a;— >7r 

whereas for initial data on another edge, a^, of the corresponding condition is 

lim'0(x) = — 00 (D.3) 

The set U~ is also infinite curvilinear strip related to Upp by the symmetry with respect 
to the axis ip. 

Let V be v-strip situated in an island Di between two v-curves and . V is 
curvilinear quadrangle bounded by and (3^ and two more bounds lying on a+ and 
a~ . Taking into account (lD.2p and (ID.3P one concludes that TV is an infinite curvilinear 
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strip stretching along U~ having the edges T^^ and T/3r. This means that TV crosses 
all a^, j = 1, . . . , at least once and pass through all the sets Dj, j = 1, . . . , N. 

Let a pair («, j) be fixed. Assume that the curves /3f are graphs of monotone non- 
decreasing functions. Then are also graphs of monotone non- decreasing functions. 
Let for all p G T~^Dj n Di the conditions gi{p) > and (72 (p) > hold. It follows from 
Lemma mi that for all the points p G T~^Dj fl Di only one of the conditions [Al]-[A4] 
hold. This means that the images TjSf fl Dj and T/3~ fl Dj consist of one connected 
component. In fact, if TjS^ crosses or a J twice one can choose two pairs of points 

such that their images qu = Tpk G Dj, 

qi = (01, q2 = (02, ^2), ^3 = (03, 03)' ?4 = (04, (p'i) 

are mismatched in the sense that at least one of the products 

(02 - 0l) (04 - 03) or (02 - 0'i) (04 - 03) 

is negative. By Lemma HI 1 the images TI3f fl Dj are graphs of non- decreasing or non- 
increasing 7-Lipschitz functions. Since Dj is an island, the boundaries are graphs 
of monotonic functions and by geometric reasons the monotonicity properties (non- 
increasing or non-decreasing) of /Sf" and T/Sf fl Dj are the same. Therefore T/3f n Dj 
are v-curves. They bounded the set TV fl Dj, therefore TV fl Dj is a v-strip. 

If the curves I3f are graphs of monotone non- increasing functions and for all 
p G T~^Dj n Di the conditions gi{p) < and (72 (p) < hold the proof repeats the 
reasoning given above making use of the conditions [B1]-[B4]. Theorem |Hl is proved. 
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